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INTRODUCTION 


ORIGWAL PAGE 
OF POOR QUALITY 



In the development of modern turbomachinery, problems of flutter instabilities 
and excessive forced response of a cascade of blades that were encountered have 
often turned out to be extremely difficult to eliminate (refs. 1,2). The study of 4 

these instabilities and the forced response is complicated, by the presence of mis- • 

tuning; that is, small differences among the individual blades. 

The theory of mistuned cascade behavior (refs. 3-8) shows that mistuning can ^ 

have a beneficial effect on the stability of the rotor. This beneficial effect is ^ 

produced by the coupling between the more stable and less stable flutter modes 
introduced by mistuning (ref. 9). The effect of mistuning on the forced response 
can be either beneficial or adverse. Kaza and Kielb (refs. 5-8) have studied the 
effects of two types of mistuning on the flutter and forced response: alternate 

mistuning where alternate blades are identical and random mistuning. 

The objective of the present paper is to investigate other patterns of mistun- 
ing which maximize the beneficial effects on the flutter and forced response of the 
cascade. Numerical optimization techniques are employed to obtain optimal mistun- 
ing patterns. The optimization program seeks to minimize the amount of mistuning 
required to satisfy constraints on flutter speed and forced response. 
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GEOMETRY OF A TUNED CASCADE 




As shovm in the figure, the blades are modeled as an Infinite cascade of air- 
foils in a uniform upstream flow with a velocity V where ^ is the stagger angle. 

* Only two degrees of freedom (bending and torsion) are considered for each blade. 

» For the tuned cascade, the blades are assumed to be in harmonic motion with h 

ar 

being the bending amplitude, the torsional amplitude and 3 the phase angle 

« between adjacent blades. For an N— blade cascade, that phase angle can take only 

• N discrete values 3 = 2irr/N. 

r 
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STRUCTURAL AND AERODYNAMIC MODEL OF BLADE 


The figure illustrates the structural model of the s-th blade. All the length 
quantities are non-dimensional with respect to the semichord, b. The blade bending 
and torsional stiffnesses are modeled through the springs and respectively. 

s s 

The effects of centrifugal stiffening due to rotation of the disk are included in 
the spring constants. The elastic coupling between bending and torsion due to pre- 
twist, shrouds and rotation is modeled through the offset distance (x^ ) of the 

s 

center of gravity from elastic axis which is located at a distance ba from mid- 
chord. The chordwise motion of the airfoil is neglected. 

The aerodynamic loads on the blade are the lift and moment per unit span L^ 

and due to motion and the lift and moment per unit span L^ and M^ due to wakes. 

These aerodynamic loads are calculated using Whitehead's extension of Theodorsen's 
isolated airfoil theory in the incompressible unsteady flow to account for cascade 
effects (ref. 10). 


z 
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EQUATIONS OF MOTION 


While the motion of a tuned cascade is simple harmonic with a fixed interblade 
♦ phase angle, the motion of a mistuned cascade is assumed to be a linear combination 

f of the tuned cascade modes. The equation of motion may be written as Eq. (1) where 

{Y} is a vector of complex amplitudes of the tuned cascade modes, [e] is the modal 
matrix containing all the possible inter-blade phase angle modes, {Q} is a forcing 
^ vector that depends on the aerodynamic wake forces and is a reference frequency. 

«■ 

When no external loads are applied the eigenvalues of the matrix [p] are cal- 
culated for a range of reduced frequencies and the flutter speed is found from the 
condition that the real part of the eigenvalue is zero. For the forced response 
calculation, the frequency and mode of the excitation has to be. assumed. In the 
present work an entire range of frequencies is scanned for the most critical forced 
response. The mode of excitation has its only non-zero component in the (N-l)th 
harmonic. 


[[P] - [I] y] {Y} = -[£]■'' {Q} 


Stabi 1 i ty 

[[P] - [I] y] {Y} = {0} 


-2N X 2N Eigenvalue problem 


0 ) 


( 2 ) 



+ i v-i 


Forced Response 

{X} = -[E] [[P] - Tj [!]]■' {Q) 


(3) 
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NONDIMENS lONAL FLUTTER SPEED, Vc/bw 


EFFECTS OF MISTUNING 


The figure illustrates the effect of alternate blade mistuning on the flutter 
speed. The flutter speed increases monotonically with an increase in alternate 
blade mistuning level. Alternate blade mistuning can have either a beneficial or 
adverse effect on forced response, depending on the harmonic of excitation. 



DESIGN FORMULATION 


The objective of the present study is to minimize the amount of mistuning 
required to satisfy given constraints on the stability and forced response of a 
cascade. The design variables are the amounts of mistuning in the individual 
blades e^, and the objective function is the sum of the individual mistunings 

raised to some integer power p. A high value of p corresponds to minimizing the 

maximal blade mistuning while p=2 corresponds to minimizing the root mean square 

mistuning. The flutter constraint is based on the result of ref. 9 which showed 

that maximum stability is obtained when all eigenvalues have the same real part y. 

Therefore, the flutter constraint is a limit on the amount of spread of the real 

part of the eigenvalues about their average value, y , as well as a requirement 

av 

that all real parts are stable. 


Design Variables: 


= Y 


al 




0 , 1 , 


N-1 


Objective Function: 


n-1 p 
F({e}) = E e- 

i=0 


Flutter Constraint: 


Spread Constraint: 



a > 0 

0 — 


Stability Constraint: 



0 


Forced Response Constraint: 


1 - -! — L > 0, 


max 


“j ^ ^c 


i = 0, 1, ... N-1 
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EFFICIENT FORCED RESPONSE CONSTRAINT 


The constraint on the forced response (Eq. (4) in the figure) requires the 
calculation of the forced response for a range of frequencies. One way of checking 
whether any violation occurs in the required range is to evaluate the response at a 
grid of frequencies dense enough to preclude the possibility of substantial con- 
straint violations between grid points. From the standpoint of the optimization 
procedure this is very costly because a constraint on the response has to be applied 
at each grid frequency. For the cases reported here a grid of 101 frequencies had 
to be used in the range 0.95 ^ 1*05. 

Two alternative techniques were used to reduce the cost of calculating the 
response constraints and their derivatives. The first is identifying local peaks of 
the response (as a function of o)) and enforcing the constraints only at these 
peaks. The main savings of this technique is in terms of derivative calculations. 
The second technique is based on the assumption that the response is most critical 
at the eigenf requencies of the stability problem. This technique results in even 
larger savings. The results in terms of number of constraints and computer time 
for a full optimization are shown in the table for a seven-blade cascade. 


RESPONSE AT 0.5 PERCENT ALTERNATE DESIGN (K=0.66) 



FREOUENCT 


LEGEND: I -*—+—+ BLADE 1 
BLADE 5 


BLADE 2 
BLADE 6 


BLADE 3 BLADE 4 

BLADE 7 


Forced Response Constraint 


1 




max 


0 
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NEWSUMT OPTIMIZER 


The optimization program used to obtain numerical results is the NEWSUMT 
program (ref. 11)* It uses the sequence of Unconstrained Minimization Technique 
(SUMT) with an extended interior penalty function (ref. 12) to represent the con- 
straints. Each unconstrained minimization is performed by using Newton’s method 
with approximate derivatives (ref. 12). The optimization procedure is particularly 
efficient when the complexity of a problem is in the constraint and the objective 
function is fairly simple. For this reason the amount of mistuning is optimized 
subject to a constraint on the response, rather than optimizing the response subject 
to a constraint on the amount of mistuning. 


Minimize F(t) 


subject to g. (t) >0, j 

J 


2, ...» m 


• SUMT - Sequential Unconstrained Minimization 

Technique 

. Extended Interior Penalty Function 

• Newton's Method with Approximate Second 
Derivatives 


OPTIMIZATION SUBJECT TO FLUTTER CONSTRAINT 


The first set of results were obtained for a twelve blade cascade subject to a 
flutter constraint of the form 

^cr 

— - 0.584 > 0 
^av 

where u is the real part of the least stable eigenvalue and y the average real 
cr 

part. The results are summarized in the table. They show that the optimized mis- 
tuning pattern is about 78.52% better than the alternate mistuning design which 
satisfies the same constraint. The maximum individual blade mistuning is 0.91% 
versus 1.4% for the alternate mistuning, and the optimized pattern is still alter- 
nating in sign. 


TABLE 1 

Results of Optimization with Flutter Constraint (k=0.66) 



Alternate Mis- 
tuning pattern 

Optimized Mis- 
tuning pattern 

Obj active Function 

23.52 X 10 "^ 

5.053 X 10"^ 

Max . mistuning 

1 . 4000 

0.9097 

E (percent) 

max'^ ' 



Least stable 

-0.002525 

-0.002526 

eigenvalue 

Mistuning 

(percentage) 


0.7628 

^1 

1.4000 

^2 

-1.4000 

-0-4768 


1.4000 

0.9018 

0 

-1.4000 

-0.6401 

^5 

1 . 4000 

0.9097 

"6 
^ 7 

-1.4000 

-0.4683 

1.4000 

0.6619 

^8 

-1.4000 

-0.6708 

^9 

1.4000 

0.2620 

£ 1 ^ 

-1.4000 

-0.7796 

10 



®ii 

£ -I 0 

1.4000 

0.1575 

-1.4000 

-0.6201 

12 
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EFFECT OF OBJECTIVE FUNCTION FOKM 


The use of the sum of the squares of the individual blade mistunings as the 
objective function is equivalent to minimizing the root mean square of the mis tuning 
pattern. Another possible objective function is the maximum individual blade mis- 
tuning. This objective function has the disadvantage of having discontinuous de- 
rivatives with respect to the design variables, e^. To avoid this problem the max- 
imum-individual-blade objective function can be approximated by the sum of a high 
power of the individual mistunings. To check whether the optimized design is 
sensitive to the objective function the optimization was repeated with the sum of 
the sixth powers of the being the objective function. The results are summarized 

in the table and show the effect to be minimal for this case. 


TABLE 2 

Effects of change in Objective Function (k=0.66) 



Optimum I 
(Exponent = 6) 

Optimum I I 
(Exponent = 2) 

Objective Function 

1.698 X lO"'^ 

5.053 X 10"'^ 

Max. mistuning 

0.8993 

0.9097 




Least stable 

-0.002525 

-0.002526 

eigenvalue 

Mistuning 

(percentage) 




0.7639 

0.7628 

^2 

-0.5662 

-0.4768 

^3 

0.8927 

0.9018 

"4 

-0.5910 

-0.6401 

"5 

0.8993 

0.9097 

^6 

-0.5352 

-0.4683 

= 7 

0.6839 

0.6619 

=8 

-0.7201 

-0.6708 

= 9 

0.2476 

0.2620 

"lO 

-0.6621 

-0.7796 


r\ ^ '7PiQ 

n T C T c 



W . J. ^ / w/ 


-0.5837 

-0.6201 

12 





OPTIMIZATION SUBJECT TO FORCED RESPONSE CONSTRAINT 


The optimization was repeated with a forced response constraint. The con- 
straint stipulated that the maximum response amplitude of the optimized design does 
not exceed the forced response of a 0.5% mistuning alternate-mistuning design. The 
two designs are compared in the table. It is shown that the objective function was 
reduced by 70% which corresponds to 45% reduction in the root mean square of the 
mistuning. 

An attempt to obtain an optimal design subject to both flutter and forced 
response constraint indicated that the alternate mistuning design cannot be improved 
upon in both categories. That is, improvements in stability resulted in deterior- 
ation in forced response, and vice versa. 


TABLE 3 

Results of Optimization with Forced Response Constraint 

(k=0.8) 



Initial Design 

Optimized Design 

Objective Function 

3.0 X 10"'^ 

9.144 X 10 '^ 

Max. mistuning 

0.5 

0.4804 

£ (percent) 



Least stable 

-0.00138 

-0.00003 

eigenvalue 



Mistuning 



(percentage) 




0.5 

0.2457 

*2 

-0.5 

-0.4804 

"3 

0.5 

0.2677 


-0.5 

-0.3049 


0.5 

0.2299 

^6 

-0.5 

-0.2620 

"7 

0.5 

0.2463 


-0.5 

-0.2024 

"9 

0.5 

0.2419 

®10 

-0.5 

-0.1170 

"11 

0.5 

0.2474 

"12 

-0.5 

-0.1123 
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CONCLUDING REMARKS 


An optimization procedure for finding optimal mistuning patterns for cascades 
subject to flutter and forced response constraints has been developed. The pro- 
cedure is based on an extended interior penalty function algorithm and seeks to min- 
imize the amount of mistuning required to satisfy the constraints. An efficient 
form of the forced response constraint which reduces computation costs by an order 
of magnitude has also been developed. 

The optimization procedure has been applied to the design of a 12-blade 
cascade and the resulting designs compared to alternate mistuning designs. It was 
found that mistuning amplitudes could be substantially reduced without hurting 
either the flutter or the forced response characteristics. However, it was not 
possible for the example problem to Improve on the alternate design subject to both 
constraints. 

The designs obtained by the optimization procedure are not practical because 
they require many different blades. Work is under way to obtain similar results 
with only 3 or 4 different blade types. 


• Optimization Procedure for Design Under 
Flutter & Forced Response Constraints 
Developed 

• Efficient Forced Response Constraint 
Resulted in Large Computer Time Savings 

• Optimized Designs Superior to Alternate 
Designs if only Flutter or only Forced 
Response is Critical 

• Number of Different Blades Should be 
Reduced 
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